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Precise measurements of the deuteron vector analyzing power AjJ and the tensor analyzing power 
A yy of the 1 H(d, 7) 3 He-capture reaction have been performed at deuteron energies of 29 MeV and 
45 MeV. The data have been compared to theoretical state-of-the-art calculations available today. 
Due to the large sensitivity of polarization observables and the precision of the data light could be 
shed on small effects present in the dynamics of the reaction. 
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I. INTRODUCTION 

Experimental studies of the three-body system are 
particularly interesting as it has become possible to solve 
the Schrodinger equation of the three-body system for 
both the ground- and the continuum states. Such three- 
body calculations, based on different techniques, are 
available from several groups 0, 0, 0> S IS 01 • Data for 
d-p elastic scattering, two-, and three-body breakup, 
and radiative capture reactions (or its inverse reaction, 
photo-disintegration) now allow for a precise quantitative 
comparison with theory. 

Combined with these calculations the radiative capture 
reaction provides an especially attractive framework as 
the electromagnetic interaction is a well understood pro- 
cess. Capture reactions in few-body systems have been 
studied for quite some time and provided, even at very 
low energy pj, valuable information on the dynamics in 
these systems. In the two-body system the cross section 
data at low energy indicated the importance of mesonic 
degrees of freedom @ and the forward angle cross sec- 
tion of the two-body photodisintegration could only be 
understood when accounting for relativistic effects 

The important role of of meson exchange currents 
(MEC) can be observed in many electromagnetic observ- 
ables. The electromagnetic form factors of the A=3 sys- 
tem as measured in elastic electron scattering jlf| or the 
electrodisintegration cross section of the two body system 
|TT| can only be understood when MEC's are taken into 
account. In these observables MEC's are intimately re- 
lated to the nucleonic S-D transitions as they give effects 
of similar size, but opposite sign. Thus, the quantitative 
study of MEC-effects requires a precise knowledge of the 
nucleonic S-D transition. The effect of such transitions 
can be enhanced in measurements of polarization observ- 
ables as the S-S amplitude, which usually dominates un- 
polarized cross sections, is strongly suppressed and al- 
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lows for a study of the small amplitudes. Based on the 
following arguments one can expect that measurements 
of polarization observables in capture reactions provide 
insight into the different roles played by nucleonic and 
mesonic degrees of freedom. 

In the energy range of 10-50 MeV the 1 H(d, 7) 3 He- 
capture process is dominated by the electric dipole tran- 
sition (El). Although MEC's can give large contributions 
to the El-transition, they can be taken into account im- 
plicitly when performing calculations with operators de- 
rived using Siegert's theorem. Thus, due to the angular 
dependence of the El-El contribution (sin 2 (# c . m .)), po- 
larization observables at medium range reaction angles 
(6c.m.) are little affected by the explicit contribution of 
MEC's. On the other hand MEC's have to be calculated 
explicitly in the magnetic transitions, particularly Ml in 
which they may give contributions up to 50% 12] . Con- 
tributions of Ml and thus MEC's are particularly large at 
small and large reaction angles. At these angles the El- 
Ml interference term, in which the small Ml amplitude is 
enhanced via El becomes dominant. Thus, angular dis- 
tributions of polarization observables with sensitivities 
to both the electric- and the magnetic transitions can be 
expected to offer excellent windows on small components 
of the interaction and provide a unique experimental test 
of MEC's. 

Measurements of polarization observables for the 
1 H(d, 7) 3 He-capture reaction with adequate precision 
are rather scarce. Difficulties associated with polarized 
beam and/or target production, and with the measure- 
ment of spin observables often led to experimental uncer- 
tainties which do not permit a significant check of the- 
oretical predictions. Although the techniques are well 
under control today, one can not expect a significant in- 
crease of the data base as the required experimental fa- 
cilities are no longer at hand [l3|. A rather complete 
account of the existing data is given in the publication 
by Anklin et al. which also discusses the results of 
our previous work in this area. 

Here we report on new measurements of vector- and 
tensor analyzing powers of the d — p capture reaction in- 
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duced with a polarized deuteron beam. In the present ex- 
periment we have extended the measurements by Anklin 
et al. in three ways. One extension concerns more ex- 
treme angles at a deuteron beam energy of 45 MeV for a 
larger sensitivity to the magnetic transitions. A second 
extension concerns the measurement of a complete angu- 
lar distribution of the polarization observables at a beam 
energy of 29 MeV. This is particularly important for the 
tensor analyzing power A yy at intermediate angles as the 
beam energy dependence shows a maximum at 29 MeV 
while at 45 MeV it is close to zero (see figure ^| . As 
will be discussed in the last section the tensor analyzing 
power at these two energies has rather different sensitiv- 
ities to the underlying physics. As a third extension we 
also measured the data of the deuteron vector analyzing 
power A% at the same kinematical points. 



to minimize systematic errors. 
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FIG. 1: Beam energy dependence of A™at 8 c . m . = 90deg. 
The data points are from Hi H dill [3 El- The dashed 
line represents the Faddeev calculation by Golak et al. 



II. EXPERIMENT 

The experiment was performed at the Philips injector 
cyclotron of the Paul Scherrer Institute (PSI) in Villigen 
(Switzerland). It used a polarized deuteron beam which 
was prepared in the PSI atomic beam ion source |19| . 
The source was equipped with a 30K cold atomic beam 
dissociator, two sextupole fields to focus (defocus) the 
atomic electrons from 2 H with spin up (down), a set of 
two strong (SF1, SF2) and one weak field radio frequency 
(RF) -transition units to induce the nuclear polarization, 
and an electron-cyclotron resonance (ECR) ionizer. The 
RF-units selected transitions between different Zeemann 
levels of the 2 H-atom hyperfine structure in an external 
magnetic field. Depending on the combination of active 
RF-units the nominal nuclear vector and tensor polar- 
ization as listed in table |U could be prepared. In the 
present experiment all five modes have been used, the 
source being cycled through them with a rate of 0.3 Hz 
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-1 
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TABLE I: Polarized beam source modes and nominal polar- 
ization values. 

Following the source, the beam was deflected into the 
injector cyclotron, accelerated to energies of 29 MeV and 
45 MeV, respectively, and guided to the experimental 
area. A schematic view of the experimental setup in the 
area is shown in Figure El 

The beam passed first a scattering chamber, where 
deuterons scattered elastically from a thin carbon foil. 
The distribution of their time of flight relative to the 
radio-frequency (RF) signal from the cyclotron was mea- 
sured using a fast plastic scintillation detector placed at 
30° below the beam axis. The cyclotron was tuned for a 
minimal beam burst width of typically 1.5 ns FWHM. 
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FIG. 2: Schematic overview of the beam line in the experi- 
mental hall NE-C. Here C(d,d) - carbon scattering chamber, 
Q - quadmpole doublet, He-pol - 4 He-polarimeter, TC - target 
chamber, D - dipole magnet to separate 3 He and unscattered 
deuterons, FC - Faraday cup. 

In a second target chamber (He-pol), a polarimctcr 
was employed for measurements of the polarization of 
the beam at regular intervals. Details of the polarization 
determination will be explained in section III Al Finally 
the beam was refocused onto a liquid hydrogen target 
mounted in the third scattering chamber (TC) for the 
measurement of the capture reaction. Downstream of 
the main target station a C-shaped dipole magnet (D) 
separated the beam from recoil particles of the capture 
reaction. A Faraday cup (FC) stopped the beam and 
measured the current. 



A. Polarimetry 

Elastic 4 He(d, a)-scattering was used to measure the 
absolute deuteron beam polarization. In the energy range 
of this experiment at a center of mass angle (i? c .m.) of 150 
deg the analyzing powers are high and precisely known 
|2Cj . In table [H] the specific values used here for the 
polarization determination of the beam are listed. 
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FIG. 3: Energy spectra at $ a i a b = 15°, Ed =29 MeV (left) and Ed =45 MeV (right). The integration limits are shown as 
vertical lines. 



Beam Energy 




A 

™-yy 


29 MeV 


0.846 ± 0.020 


0.910 ±0.016 


45 MeV 


0.497 ±0.011 


0.921 ±0.013 



TABLE II: Vector (A a ) and tensor (A yy ) analyzing powers 
for elastic d — a scattering at a center of mass angle of 150 
deg. 



A 4 He gas cell with 5 /zm thin Havar windows, oper- 
ating at a pressure of 50 kPa, served as the target of the 
polarimeter. At laboratory angles of ±15° correspond- 
ing to t? c . m . = ±150deg recoil a-particles were detected 
with two symmetrically arranged "passivated implanted 
planar silicon" (PIPS) detectors. A double slit system 
with Ta-collimators at distances of 12 cm and 45 cm from 
the center of the scattering chamber shielded background 
particles from beam-target window reactions. The colli- 
mators limited the angular acceptance to ±0.5 deg. PIPS 
detectors with thicknesses of 700 /im (500 /mi) at beam 
energies of 45 MeV (29 MeV) allowed for the discrimina- 
tion between the stopped recoil a-particles and the high 
energy elastically scattered deuterons which deposited 
less energy as compared to the stopped a-particles. Fig- 
ure[3]displays two sample energy spectra measured at the 
two deuteron energies. 

The polarization has been determined with the expres- 
sion for the cross section of a reaction with a mixed- 
polarization deuteron beam: 

(s)^(s)„( 1 + l^ + ^») w 

with i = b,c,d,e. Vector (tensor) polarizations of the 
beam for source state i = 6, c, d, e are denoted with p\ 
(Pzz) ( see table HJ and (da/dYl) a is the unpolarized cross 



section. Coordinate system and symbol definitions fol- 
low the Madison convention 21]. Based on the general 
expression the following equations can then be written 
for each polarization state i: 

,KiNi-Nl KiNi-N°, 1 



4 _ KjNj-Nl KiNL-Nl 1 
Pz 1 N° N° >3A y {) 

with Nl. (NL) the number of events detected by the left 
(right) detector and K\ [K % _) correction factors which 
account for dead time (DT) and Faraday cup (FC) dif- 
ferences for the different states i. Cross sections can be 
replaced by accumulated counts because solid angles and 
efficiencies cancel. 

The polarization of the deuteron beam was measured 
every 8 hours. As already observed in the previous ex- 
periment by Anklin et al. |14| the polarizations of the 
source have been very constant over a run period of 14 
days. Typical mean polarization values were p\ = 0.25 
for the vector and p l zz — 0.65 for the tensor polarizations 
for states i = b — d with accuracies of 2-3%. 

B. Setup for the 1 H(d, 7) 3 He analyzing power 
measurements 

Figure^lgives a detailed view of the setup for the mea- 
surements of the 1 H(d, 7) 3 He-capture reaction. The po- 
larized deuteron beam was incident on a liquid hydrogen 
target (LH-T) with a thickness of 14 mg/cm 2 enclosed 
by 2.5 fim Havar windows. The target cell was cooled 
to about 16 K with a closed-cycle helium-refrigerator 
and operated at a pressure of 0.25 bar. It was mounted 
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in a specially designed scattering chamber with 3 mm 
Al-walls and a conical entrance beam tube to allow for 
measurements at very large scattering angles. 

The capture photons were detected by four large 
BaF2-counters in coincidence with the recoil 3 He parti- 
cles detected in thin plastic scintillators (R). To separate 
the beam from the recoil particles a C-shaped dipole- 
magnet (D) with a vertical pole-tip distance of 90 mm 
was used downstream of the interaction point. 




FIG. 4: (Color online) Schematic overview of the setup: LH- 
T - liquid hydrogen target, D - dipole magnet to separate 3 He 
and unscattered deuterons, R - recoil-detectors, FC - Faraday 
cup. 



1. "/-detectors 

The photons from the capture reaction with energies 
from 13 MeV to 16 MeV for E d =29 MeV and from 
17 MeV to 24 MeV for E d =45 MeV were detected with 
four BaF2-scintillators. The detectors were placed at a 
distance of 80 cm from the target at various angles in the 
range from 27° to 169°. Each detector consists of four 
large cubic crystals 8x8x25 cm 3 placed in aluminium 
containers with a wall-thickness of 5 mm. The boxes 
were shielded with 5 cm of lead on the sides and 5 cm 
borated plastic in front. 

The scintillator material was chosen because of its ex- 
cellent timing characteristics. The light response of BaF2 
is characterized by two decay times of 0.7 ns (short) 
and 620 ns (long) and a light emission spectrum with 
maxima at 220 nm and 310 nm, respectively [2^|. The 
BaF2 fast component allows for an excellent discrimina- 
tion between the capture photons and the copious neu- 
trons from break-up reactions and beam-target window 
interactions. This is particularly important for the mea- 
surements at extreme angles where the capture cross sec- 
tion is very small. Due to its high density of 4.89 g/cm 3 
BaF2 has a very good efficiency and an acceptable en- 
ergy resolution of about 16% (short component) for E^xs 
20 MeV. Each crystal is connected to fast photomultiplier 
tubes Philips XP4318B with special optical gel trans- 
parent for UV-radiation (Baysilone Ol M 600000 by GE 



Bayer Silicons GmbH). A light emitting diode, operated 
at a frequency of 10 Hz, allowed to monitor the gain 
during the experiment. At forward angles the standard 
/i-metal shielding of the photomultiplier tubes has been 
complemented with two closed /i-metal boxes with a wall- 
thickness of 1.5 mm, shielding the whole detectors from 
the magnetic field of the C-magnet of about 50 G. 



2. Recoil- detectors 

To detect the deflected recoil particles with energies 
of 17-21 (26-31) MeV at the beam energy of 29 MeV 
(45 MeV) plates of plastic scintillator Pilot U with a 
thickness of 1.2 mm were used. The 3 He recoiled in a 
cone between 0.4° and 2.6° at both deuteron energies de- 
pending on the photon angle. The strength of the mag- 
netic field was chosen to deflect the initial deuteron beam 
by ~ 8° such as to separate the 3 He and deuterons by 
more than 10°. 

The position and size of the recoil detectors were de- 
fined by numerical simulations of deuteron and 3 He tra- 
jectories in a magnetic field based on measurements of 
the magnetic field map determined at PSI. Angular and 
spatial deviations due to beam divergence, beam size and 
multiple scattering in the target was taken into account. 
The detectors were designed to accept 99% of the recoil 
particles. The size of the individual Pilot [/-pieces was 
chosen to yield nearly equal 3 He fluxes in all detectors to 
protect the electronics from signal over-load. The thick- 
ness was chosen by the requirement to stop the recoil 3 He 
but not the Rutherford scattered deuterons and the pro- 
tons from break-up reactions. The photomultiplicrs were 
mounted vertically at a distance of about 50 cm where 
the magnetic field from the deflection magnet could be 
sufficiently shielded. 

3. Electronics 

All the detector signals were multiplexed to form a 
trigger signal and to process the signal for time and am- 
plitude measurements. The trigger signals were clipped 
to correct for base line shifts at high rates and fed to 
constant fraction discriminators to minimize the signal 
amplitude dependence. A coincidence was requested be- 
tween each recoil and BaF2 detector within a time win- 
dow of 25 ns. The sum of all coincident signals was used 
for a further coincidence with the RF-signal of the cy- 
clotron. To form these coincidences first is important to 
minimize electronic dead time effects at the high rates. 
The final trigger signal with the RF time provided the 
start of the TDC and triggered the read-out of the CA- 
MAC system. For each detector the retimed coincidence 
signal was used for the gate of the charge integrating 
ADC. For redundancy coincidences between the recoil de- 
tectors and the BaF2 detector were also formed. For the 
BaF2 detectors two amplitudes were recorded. One with 
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a short 25 ns gate for the fast component and one with a 
700 ns gate for the long one. For a time of flight (TOF) 
measurement, the TDC's were started with the beam RF 
signal and stopped with the individual fast component 
detector signal. The same was done for the recoil de- 
tector signals. A signal from the digitized beam current 
of the Faraday cup, a real-time clock and pulser signals 
as well as the current polarization state signal were fed 
into scalers. These signals were relevant to correct false 
asymmetries from dead time and beam current variations 
correlated with the polarization state. All the informa- 
tion was written on an event -by-event basis on disc for 
on-line analysis and also written to tape for backup and 
replay. 

III. ANALYSIS AND RESULTS 



software time cuts on the prompt 7 events. In order to 
achieve this result software cuts have been applied on 
the BaF 2 -recoil-TOF, the BaF 2 -RF-TOF, and on the 
energy information of the recoil-detectors. In addition, 
due to the kinematical angular correlation the proper re- 
coil detector is selected for a given BaF 2 . This selection 
provides for a significant reduction of background. 

Figure shows an example of a two-dimensional his- 
togram with a BaF 2 detector lightoutput versus its 
BaF 2 -recoil-TOF of the corresponding recoil detector. 
The cuts on the BaF 2 -RF-TOF and recoil lightoutput 
are applied. The two-dimensional histogram shows a 
clear separation between the 7— 3 He coincidences and the 
unstructured accidental ones. With an additional cut on 
the coincidence time BaF 2 -energy spectra like the one 
shown in figure El have been achieved. 
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FIG. 6: (Color online) Example of a 7-energy spectrum with 
all cuts applied and the integration limits for the capture 
events. For studies of the background analyzing power as 
discussed in the text, a third integration line including the 
background is also shown 



FIG. 5: Example of a spectrum of the short component am- 
plitude of the BaF2 response versus the TOF between the 
7-rays and the recoil particles. Cuts are applied on the recoil 
energy and on the BaF 2 -RF-TOF. 

The main challenge of the data analysis is to single 
out the few capture events from a huge background due 
to hadronic reactions. Particularly for data taken at the 
very large and very small scattering angles the applied 
hardware coincidence is indispensable. The coincidence 
requirement results in a signal to noise ratio of the or- 
der of 1 for these data. The remaining background is 
mostly due to accidental coincidences and n — p breakup 
reactions of deuterons on hydrogen and the target win- 
dows. Depending on the kinematics the signal to noise 
ratio could be enhanced to at least 25 due to the excellent 
timing resolution of the BaF 2 detectors and the efficient 



In figure El a 7-energy spectrum (BaF 2 short compo- 
nent) with all mentioned cuts applied (except the cut on 
BaF 2 -lightoutput) is shown. The spectrum also shows 
the integration limits applied to determine the capture 
events. The left part of the spectrum is due to the re- 
maining low energy 7 background which also contributes 
to the region of the 7-peak within the integration lim- 
its. To determine the effect of this remaining background 
contribution the response function of the BaF 2 -detector 
to monoenergetic photons must be known. 

For this determination an additional experiment with 
monoenergetic 7-rays was performed at the Physics Insti- 
tute of the University of Basel. Monoenergetic 20 MeV 
7-rays were produced in a 3 H(p,7) 4 He reaction with a 
1 MeV proton beam provided by a Cockroft- Walton ac- 
celerator. The detector was placed at 110 deg close to 
the maximum of the angular distribution of the photons. 
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FIG. 7: (Color onlie) On the left the short component of the BaF2 response function for 20 MeV 7-rays. On the right 
experimental 7-spectrum in comparison with the fit-function. The integration limits are represented by two vertical lines. The 
contribution of the exponential function and the folded peak function are also shown separately 



The distance between the tritium target and the BaF2- 
detector was 80 cm in order to reproduce the geome- 
try of the d — p-capture experiment. The target was a 
0.45 mg/cm 2 Ti layer on a Cu-plate with 2 Ci absorbed 
tritium. The electronics was a simplified version of the 
d — p-capture experiment using the same electronic mod- 
ules and gate widths. The measured short component of 
the BaF2 response function for the monoenergetic pho- 
tons from this study is shown on the left in figure 

The sum of the experimentally determined response 
function and an exponential function folded with a Gaus- 
sian is used in the analysis to fit the spectra of the short 
component of the BaF2 detectors. The peak position, the 
amplitude, the width of the Gaussian and the slope of the 
exponential function are used as parameters. An exam- 
ple of such an analysis is shown on the right in figure 
The reduced \ 2 °f the fits varies from 0.8 to 2.2. The 
upper end of the peak is distorted due to pile-up during 
the d — p-capture experiment and thus is excluded in the 
fit. 

As the energy of the photons in the d — p-capture ex- 
periment varied between 13 MeV and 24 MeV the energy 
dependence of the response function was studied with a 
GEANT-Monte-Carlo simulation. Photons with energies 
of 13 MeV or with 24 MeV incident on a BaF2 crystal 
folded with a Gaussian to fit the experimental peak width 
have been compared to a simulation at 20 MeV, the en- 
ergy of the model peak. Scaling the response functions 
to the model peak energy resulted in an energy depen- 
dence of the response of less than 5% for the integration 
limits applied in the analysis. This results in a negligi- 
ble relative error contribution of 0.2% for a background 
contribution of 4.6%. 

With the assumption of an unpolarized background the 
contribution from 0% to 4.6% within the integration lim- 



its can be interpreted as a dilution factor N tot / (N tot — 
jybackgry rp Q accoun ^ f or this dilution the corrected ana- 
lyzing powers are then given as 



a capture a extracted 

v, yy ~ v, yy 



J\J~tot jybackgr 



(3) 



Here N tot is the total number of events within the inte- 
gration limits and N backgr is the total number of counts 
of the fitted exponential function within the integration 
limits. Ay Xt yy Cted are the vector-, tensor-analyzing pow- 
ers determined from the total number of counts in the 
different polarization states, respectively. 

In order to verify the assumption of an unpolarized 
background the following tests have been performed. 
First, the corrected analyzing powers at a given kine- 
matic point must be independent of the dilution factor. 
Different sets of runs with different background condi- 
tions due to different beam currents have been compared. 
Within the statistical errors no systematic deviations 
have been found. Second, if the analyzing powers of the 
background would not be zero, they would contribute to 
the extracted analyzing power values as: 



A capture 

y, yy ~ 



ATtot 

A extracted _ 

V, yy ' jytot _ jybackgr 



\ back 



jyba 



back 



Vi VV jytot _ ]\Jbackgr 



(4) 



Large background contributions with finite analyzing 
powers would significantly distort the corrected capture 
values. This can be checked by artificially increasing in- 
tegration limits in order to include a large background 
fraction in the low-energy part of the spectrum (see figure 
EJ. The corrected analyzing powers from these integra- 
tion limits can be compared to analyzing powers within 
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integration limits with essentially no background contri- 
bution. It is found that the analyzing powers calculated 
with the larger integration limits are completely consis- 
tent with the analyzing powers calculated for the capture 
peak. This confirms that within the statistical limits the 
background contribution is unpolarized, hereby justify- 
ing the correction described with equation [3J 

The resulting deuteron vector (A^) and tensor (A^,) 
analyzing powers from the capture reaction are given in 
table lllll and llVI The statistical and the systematic errors 
are listed separately. As can be seen the statistical errors 
dominate the total error for all data points measured. 



u c.m. (d — 7) 


A y ± SA y tat ± 5A" y ya 


A 4- A A atat + A A stat 


deg 


xlOO 


xlOO 



33.52 -3.821 ± 0.335 ± 0.084 3.010 ± 0.383 ± 0.090 

55.32 -1.808 ±0.252 ± 0.061 2.197 ± 0.288 ± 0.070 

76.41 -0.443 ± 0.262 ± 0.050 2.236 ± 0.300 ± 0.053 

116.16 0.465 ±0.253 ±0.051 2.452 ± 0.289 ± 0.054 

134.94 1.030 ± 0.285 ± 0.052 2.182 ± 0.326 ± 0.055 

153.19 0.325 ± 0.470 ± 0.044 0.438 ± 0.538 ± 0.045 

170.21 1.097 ±0.412 ±0.060 -2.934 ± 0.520 ± 0.067 



TABLE III: Vector (A y ) and tensor (A yy ) analyzing powers 
for the d — p radiative capture reaction at a deuteron energy 
of 29 MeV. 



u c.m. (d — 7) 


A y ± SAf at ± SA y ya 


A 4- A A stat + A A stat 
^-yy yy yy 


deg 


xlOO 


xlOO 



31.04 -5.134 ±0.186 ±0.095 

51.20 -1.955 ±0.160 ±0.042 

112.88 1.837 ±0.190 ±0.042 

137.89 2.986 ±0.139 ±0.057 
170.48 4.507 ± 0.328 ± 0.081 



4.610 ±0.228 ±0.113 
1.199 ±0.197 ±0.048 
1.891 ±0.236 ±0.049 
2.014 ±0.171 ±0.068 
1.334 ± 0.407 ± 0.095 



TABLE IV: Vector (A y ) and tensor (A yy ) analyzing powers 
for the d — p radiative capture reaction at a deuteron energy 
of 45 MeV. 



Various systematic errors are accounted for. A relative 
10% error is estimated for the background contribution. 
In addition, correlations between luminosities and polar- 
ization states could lead to false asymmetries. To deter- 
mine a systematic error the asymmetries of two unpolar- 
ized runs with different current have been used to deter- 
mine a false asymmetry due to different dead times and 
luminosities. The determined effect can be scaled to the 
differences in currents and dead times for different polar- 
ization states present during normal running conditions. 
A resulting systematic uncertainty of 0.00043 (0.00024) 
results due to current differences of 0.005(0.003) nA dur- 
ing the data taking at _Ed=29(45) MeV. Including the 
uncertainty of the polarization determination the total 
systematic errors vary between 0.00044 and 0.00113 com- 
pared to the statistical errors from 0.00187 to 0.00470. 



IV. COMPARISON TO THEORY 

The present data together with previous data taken at 
the same deuteron energies are c omp ared to three dif- 
ferent recent calculations [23L l2~il |25| . The calculations 
are all exact in the sense that they provide a full solu- 
tion of the Schrodinger equation from a realistic nucleon- 
nucleon interaction for both the ground- and continuum 
states. An exact treatment using the same Hamiltonian 
for initial and final state is known to be essential for a 
successful description of polarization observables of the 
d—p capture reaction. Whereas for the description of the 
differential cross section the use of the Born approxima- 
tion leads to an underestimate of order 15%, the tensor 
analyzing power A yy is underpredicted by more than a 
factor of 2. The importance of initial state interaction in 
these observables has been observed already in the cal- 
culation by J. Torre More recently a calculation 
by A.C. Fonseca and D.R. Lehman |2(| confirmed that 
re-scattering effects in the initial state are crucial for a 
determination of polarization observables. 

The techniques applied in the calculations by the three 
gro ups are very different. The Bochum-Cracow group 
[23j obtains the wave functions by solving the Faddeev 
equations in a non-relativistic framework. In the cal- 
culation of the PISA group |2J| the wave functions are 
calculated using the pair-correlated hyperspherical har- 
monics method |4|. Both groups use the Argonne vl8 
two-body potential [27J as underlying nucleon-nucleon 
potential. In addition, the Urbana IX three-body poten- 
tial [2^| is also included by both groups. In the approach 
by the Hannover group |25j | the three-particle scattering 
equations are calculated exactly with a Chebyshev expan- 
sion of the two-baryon transition matrix. This approach 
employs the CD-Bonn potential (2^1 as the underlying 
two-body-potential. The three-nucleon-force (3BF) is 
generated here via a coupled channel extension of CD 
Bonn which allows for a single nucleon transition to a 
static A isobar. The CD-Bonn + A extension is as ex- 
act as CD-Bonn as it is also fitted to the experimental 
two-nucleon data up to 350 MeV |30j . 

In figure |H1 the results of the three calculations for 
the analyzing powers at the energies of the experimental 
data are compared. In this comparison only one-body- 
currents arc included in the calculations. In addition, the 
results of one calculation by Skibinsky et al. (dash) is 
shown which neglects the 3BF. The figures confirm that 
when accounting for one-body-currents only, the results 
are essentially the same. Small deviations can be ob- 
served for the calculation which does not include a 3BF 
(dash) and for the calculation which includes the 3BF as 
a static A isobar (solid). Thus, the effects of the 3BF are 
small but notable in A^ when a different 3BF-model is 
employed. 

Figure [3] also shows the experimental results of the 
present work together with the results of and [l2| . 
Whereas the results of the three experiments are in good 
agreement the comparison to the calculations shows a 
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FIG. 8: Ay (top) and A yy (bottom) for 29 MeV (left) and 45 MeV (right) incident deuteron energy as a function of the 
center-of-mass angle between deuteron and outgoing 7. Data of the present experiment (•) together with the data by Anklin 
et al [l4T| (o) and Jourdan et al. |T3| (□) are compared to the one-body-calculations by Deltuva et al. (solid), Skibinski et al. 
with/without 3BF (dotdash/dash), and Viviani et al. (dot). 



large discrepancy. As will be shown, the dominant cause 
of these deviations are the missing two-body-currents in 
the calculations. One should note that all the theoreti- 
cal results have been folded with the acceptances of the 
experimental data of ±5.7 deg. 

In calculations of capture- and photodisintegration ob- 
servables the dominant part of many-body-currents is 
usually included using the Siegert theorem [3l[. The 
standard Siegert theorem is formulated in a multipole ex- 
pansion in which part of the transition currents can be re- 
placed by the charge operator. Alternatively is has been 
shown by Golak et al. 2] that in momentum space a mul- 
tipole expansion of the current operator is not necessary; 
this approach is applied in the calculations by Skibinski 
et al. Whereas two body currents are implicitly included 
in the dominant part of the electric transitions, they are 
not accounted for in the magnetic transitions and a small 
part of the electric ones. In the calculation by Skibinski 
et al. only the one-body-current is used in these latter 
parts of the transition. An alternative approach is used 



in which the 7T- and p-exchange currents are taken into 
account explicitly using the Riska prescription |32| . 

The effect of the two-body-currents is shown in fig- 
ure [|J] The description of the data has largely improved 
which confirms the presence of large two-body-current 
effects. The solid lines give the results including the tt- 
and /5-exchange currents explicitly whereas dashed lines 
show the results from the Siegert approach. In general, 
the description of the data is much better with the ex- 
plicit treatment of the exchange currents. This suggests 
that magnetic transitions, for which MEC's are not ac- 
counted for in the Siegert approach, are a relevant part 
of the transition. As expected, this holds particularly 
for the wings of the A yy data. In the intermediate an- 
gular range, which is dominated by the El-transition, 
one would expect the Siegert calculation, which implic- 
itly accounts for MEC's, to be the more successful ap- 
proach. However, only the A yy data at the 45 MeV are 
described well. This is surprising as these data are par- 
ticularly sensitive to small ingredients of the calculations. 
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FIG. 9: Same as figure [S] Here the data are compared to the calculations by Skibinski et al. with an explicit treatment of the 
exchange currents (solid) and within the Siegert approach (dash). 



Thus, one should conclude that only a more complete ex- 
plicit account of the two-body-currents can improve the 
description of the data. 

Such a more complete approach for many body cur- 
rents has been employed by the PISA group. This calcu- 
lation includes 7r— , /?-, w— , and a-exchange currents con- 
sistent with the NN-potential as well as the currents asso- 
ciated with the pit"/ and luttj transition mechanisms and 
with the excitation of intermediate A-resonances. One 
should note that the pir^ and lottj terms are not fitted to 
data but used with the standard coupling constants |33fl. 
In addition, these terms give very minor effects. An alter- 
native approach to the one by Skibinski et al. to derive 
the exchange currents is employed |24| . The currents are 
strictly consistent with the interaction potential. Thus, 
three-body-currents are included for calculations with 
the Urbana IX three-body potential. An essential as- 
pect of the calculation is that the current conservation 
relation is satisfied in all calculations. 

An additional improvement of the calculations by Mar- 
cucci et al is the inclusion of the point Coulomb inter- 
action. To account for the Coulomb interaction between 



the two protons is straightforward in these calculations 
as they are performed in configuration space. The ef- 
fect is small for the analyzing powers of the present work 
but, as will be discussed below it improves the compar- 
ison between data and calculation. Figure ITUI shows the 
comparison of the data with these theoretical results. 

In particular the description of the A y data for both 
energies improved with this approach in most of the an- 
gular region. If three-body-current effects are included 
the description also agrees with the A yy data at 30 MeV 
but still shows deviations at 45 MeV. In particular, the 
fall-off of the calculations in A yy at 45 MeV at backward 
angles can not be removed by many-body-currents. De- 
viations from the data are also present at very backward 
angles for the Ajj data. This suggests that additional ef- 
fects are relevant for a complete description of the data. 

A similar discrepancy was a long standing puzzle in the 
deg cross section of the two-body photodisintegration 
because relativistic effects had been considered unimpor- 
tant at these low energies. Cambi, Mosconi, and Ricci 
solved this puzzle demonstrating the importance of the 
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FIG. 10: Same as figure[S] Here the data are compared to the calculations by Marcucci et al. with 2-body currents (dash) and 
with 2- and 3-body currents (solid). The calculations by Skibinski et al. with the explicit treatment of the exchange currents 
are also shown (dot-dash). 



relativistic spin-orbit contribution to explain the discrep- 
ancies between data and calculations. 

In a (fc/mjv)-expansion of the current operator with k 
and mpf the nucleon momentum and mass, relativistic ef- 
fects of order (fc/m/y) 2 are included in the calculation by 
Dcltuva et al. |25| . The two-body currents are included 
in this calculation using the Siegert approach without the 
long-wavelength approximation often applied. In con- 
trast to the calculation by Skibinski, explicit one- and 
two-body currents in the non-Siegert parts are also ac- 
counted for. 

In addition to relativistic effects, also this calculation 
includes the Coulomb interaction between the two pro- 
tons |34| . For calculations performed in momentum space 
the Coulomb potential is usually omitted due to con- 
vergence difficulties. To include it a screened Coulomb 
potential is used, corrected for the unscreened limit us- 
ing a renormalization procedure. A recent compari- 
son with the calculations of the PISA-group discussed 
above demonstrates the reliability of the momentum 
space calculation [35| . Figure ITU shows the data in com- 
parison with the calculations by Deltuva et al. With 



this calculation the data can be reproduced over most 
of the angular range. The dotted line represents the 
results without relativistic- and Coulomb effects. The 
dashed line shows the results with relativistic corrections 
included and the results shown with the solid line include 
also the Coulomb interaction. In particular the backward 
angle fall-off of A yy can be removed with the relativistic 
spin-orbit effect. In addition, also the backward angle 
deviations of the non-relativistic calculations in A y are 
improved. 



V. CONCLUSIONS 

Precise deuteron vector analyzing powers, A^, and ten- 
sor analyzing powers, A yy , have been measured for the 
1 H(d, 7) 3 He-capture reaction at two incident deuteron 
energies. The energies have been chosen in order to em- 
phazise different dynamical effects. A yy for intermediate 
angles shows a maximum at 29 MeV, whereas it crosses 
zero around 45 MeV. 
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FIG. 11: Same as in figure[H] Comparison to the calculations by Deltuva et al. with 2-body currents (dot), plus relativistic 
corrections (dash), and plus added Coulomb corrections (solid). 



The data have been compared to modern three-body 
calculations from three different groups 0, |25| . Al- 
though the groups employ very different techniques to 
compute the wave functions for ground- and continuum 
states, the computed polarization observables agree very 
well with each other when only one-body-currents are in- 
cluded. However, large discrepancies are present between 
such calculations and the data. It has been shown, that 
the dominant part of this discrepancy can be corrected 
for accounting for two- and three-body-currents. How- 
ever, it could also be shown that relativistic effects and 
the Coulomb interaction between the protons play a role. 
The agreement between calculation and data in certain 
angular ranges improves significantly when such effects 
are included. 

In summary, the present precision data represent a 
challenging testing ground for "state-of-the art" three- 
body calculations. For a complete description of the data 



two- and three-body currents as well as relativistic- and 
Coulomb effects have to be taken into account. Due to 
their precision the data are very sensitive to various small 
effects which usually are unobservable in cross section 
data. They clearly demonstrate the power of the com- 
bination of polarization data with the electromagnetic 
process [l3j . 
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